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Abstract 

The effects of stochastic absorption and ejection of molecules by 
growing droplets have been considered. Both analytical and numer- 
ical approaches have been used. They demonstrate the satisfactory 
coincidence. It is proved that in general case corresponding to the 
asymptotic at big numbers of molecules in the critical embryo the ef- 
fects of stochastic growth are small in comparison with the effects of 
stochastic appearance of droplets. 

This paper continues the set of publications [2], [5], [7] about stochas- 
tic effects in nucleation. It is known that the growth of embryos occurs in 
nucleation kinetics stochastically. Really, in a given moment of time t the 
probability to absorb a molecule to a given cluster of a size v (i.e. contain- 
ing v molecules of a liquid phase) from the time t up to the time t + dt is 
W + dt where dt is an elementary interval and W + is the absorption coefficient 
calculated in the free molecular regime of substance exchange as 

W + = ^v tP Sa 

Here Vt is the thermal velocity of motion of the molecule of vapor, p is 
the density of vapor, S is the surface of the cluster, a is the condensation 
coefficient. This formula can be easily derived from the simple gas kinetic 
theory. 

The evident restriction for dt is dt -C (W + )~ 1 . 

There exists an analogous probability to eject a molecule into a vapor 
phase W~dt where W~ is the ejection coefficient. Ordinary one can put W~ 
equal to W + at the density of vapor equal to the saturated vapor over the 
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curved surface of the embryo. Precisely speaking this is one of the postu- 
lates of nucleation theory stating that the internal state of embryo does not 
strongly depend on the state of vapor. So, 

W~ = W + \ p=Pe(l/) 

where p e is the density of vapor in equilibrium with the embryo of the size 
v. 

Here the additional restriction for dt is dt <C (W - ) -1 . 

Since the most frequent ejections and absorptions take place for big em- 
bryos one can check restrictions for the "maximal" droplet in the system. 
Certainly, this is the supercritical droplet and here W~ < W + . So, one can 
check only dt < (W + )~ 1 . 

So, we can see that the process of growth is principally stochastic and 
required the corresponding description. 

From the first point of view to use the law of regular growth the following 
arguments are presented: 

• The characteristic size of droplets is rather big. So, the characteristic 
dispersion which is about the square root of the total sum of both 
absorption and ejection events is enough small in comparison with the 
mean number of molecules inside the cluster. 

• Ordinary the number of droplets is so giant that the averaging over the 
whole volume of the system leads to the compensation of deviations of 
particular droplets from the mean value of the size coordinate. 

Both these remarks explain the negligible character of stochastic correc- 
tions in droplets growth and we have to discuss them. 

About the first remark one can come to the following notes: 

• It is well known that the rate of embryos growth has the avalanche 
character (see pQ). Namely the avalanche consumption leads to the 
fact that the dispersion 8v of the number of molecules v inside the 
droplet doesn't satisfy the ordinary relation 8v ~ v 1 / 2 . Here, in the 
free molecular regime of growth we shall show that bv approximately 
has another power behavior bv ~ z/ 1_e with e close to 0. 

So, the relative weight of fluctuations doesn't disappear rapidly. For- 
tunately, e is positive and the formal convergence takes place. 
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• It is quite possible that due to the complex dependence of W + , W~ 
on v the mean value < v > of the droplet size v averaged over many 
attempts differs from the value v reg calculated on the base of regular 
law of droplets growth. 

So, the special consideration is necessary. Fortunately in some sense (it 
will be clear later) one can show that v Teg is close to < v >. 
The second remark faces the following notes: 

• The observed number of droplets in the system isn't giant, this specific 
feature has been fully explained in investigation of effects of stochastic 
appearance of supercritical clusters |2|. 

• Even being averaged over the giant number of droplets the mean co- 
ordinate can differ from the value predicted on the base of the regular 
rate of growth. The reason can be the non-linear dependence of u(t) 
on t. 

Here the role of only discrete values of v is also important, but even 
with account of discrete character of v one can not come to coincidence of 
approach based on the regular growth and approach of stochastic growth. 
Here the difference is very small but still it exists. 



1 The model 

The complex dependence of W~ on v occurs mainly through the dependence 
of p e on v. But one can not get precise final results in analytical form taking 
into account namely this dependence without extracting asymptotes. 

For the supercritical clusters one can use instead of the mentioned density 
the density of the vapor saturated over the plane surface. Then 

w- « w - = W+\ p=Pe[oo) 

where p e (oo) is the density of the saturated vapor over the plane surface. 

In further analysis we shall use W ~ instead of W~ . 

In renormalized scale of time which we choose for simplicity one can write 
that 

Wq = 3z/ 2/3 

w + = 3(c + iy /s 



3 



Here ( is the supersaturation of the vapor 

C = p/pipo) - 1 

The regular law growth for the supercritical cluster in this time scale can 
be written as 

^ = W+ - W ~ = 3C^ 3 

Here this law is written already for the supercritical droplets. 

Stochastic model used in numerical simulation will be the following: 

• At initial moment of time t — the droplet is situated at v — v 

• At every step dt two random numbers r + and r_ are generated. If 
r + < W + dt then v goes to v + 1, if r-i < W_dt then v goes to v — 1. 

• At t = tf in the process of growth stops and the attained value Uf in will 
be the result of calculations 

Several (many) attempts have been made and < v > is calculated as 

J2 v jin 



< V >- 



N 



where iV is the number of attempts. The value of characteristic fluctuation 
8v will be calculated as 1 



One has to stress that it is impossible to put uq close to because we 
consider the supercritical droplets. Moreover when u Q goes to one can 
expect divergence of bv. Fortunately this divergence is rather weak. It will 
be discussed later. 

Another reason to forbid the small values of u is that in this process one 
can not attain v = in principle. What shall we do with the totally dissolved 
embryo? This is the necessary disadvantage of the approach used here. So, 
one has to take z/ greater than 10. The last value ensured the necessary 
accuracy and negligible intensity of the process of dissolution. 

1 The value < v > here is calculated on the base of the separate set of attempts. 
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2 Example for the mean coordinate 



Let the characteristic size of the cluster be Uf in = 1000 molecules. Let the 
supersaturation be ( = 4 (during the nucleation period it falls but not so 
essentially, so one can take the characteristic value). 

We start at uq = 20 and try to attain Ufi n = 1000. We calculate the time 
necessary to attain it according to macroscopic regular law in continuous 
approximation 

1/3 1/3 
, _ U fin ~ V st 

tcont 



c 

(in corresponding time units) and get t cont = 1.8214. Then we recalculate the 

1 T-^i = Vfir, 

3 l^i=v at 



time according to regular law with discrete jumps (t disc = \ Y?i=jT i 2 ^ 3 



and get tdisc = 1.8267. 

The value t cont has to be corrected up to the integer number of elementary 
time intervals dt. In this example dt = 0.00002. Then Uf in will be slightly 
another. 

It is more profitable to fulfil the simple simulation of regular growth. At 
every step dt instead of v we get v + 3(v 2 / 3 dt. The final value will be marked 
as Vff. This way allows to estimate the deviation caused by the finite value 
of dt. 

The general result of simulations is that the value i/ff is practically the 
same as < v >. 

To note the difference between these values one has to see that 

| < v > -is ff \ > 5v/VN (1) 

So, we have to calculate 5u. 

Now we turn to the dispersion. The mean number of events can be 
estimated as 

< i^tot >— ^ Vf in 

So, the estimate for the mean standard deviation of the size v from the mean 
value < v > in a particular attempt will be 



5v = ^2 < E tot > 

Namely, 

Su S3 V2000 w 45 
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Numerical simulations give 5 ca iu ~ 155. 

The discrepancy will be a matter of discussions. 

For the final values we have the following results 

v f f = 1004 

< v >= 998 

To ensure that 

< V Vff ± Vf in = 1000 

we fulfilled 10000 attempts. Formally the necessary inequalities (0) are satis- 
fied. But still this discrepancy can be the result of small errors in simulation. 
Really, the random number generators are not so perfect. Also there has 
to be a special correction caused by the prescribed sequence of possible ad- 
sorption and ejection events. This weak correlation remained without our 
attention. 

Certainly, the deviation in the mean coordinate is out of practical interest. 
It is much more interesting to consider the deviation in dispersion. 

3 Analytical explanation for the mean value 

To see that < v >~ Uff we start from the master equation in Fokker-Planck 
approximation for the partial distribution function n over size v: 

at av av n e 

where n e is the formal equilibrium distribution. The last equation is valid 
for small values of supersaturation. 

From the equation of detailed balance it follows that 

W~{y + l)n> + 1) = W + {u)n e (u) 

Since 

n e ~ exp(— F(u)) 

where F is the free energy of the embryo, one can get explicit kinetic equation 
on n. 

In the cappilary approximation 

F = -ln(C+ l)z/ + a^ 2/3 
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where a is a scaled surface tension. 
For the supercritical embryos 

dF/dv « -ln(C+ 1) 

Taking into account that W + = 3(( + 1)v 2/3 /t, W~ = 3z/ 2/3 /r and acting 
in a scale where 3/r = 1 one can get 

The next step is to note that asymptotically at v — > oo 

9 ( + 2/3 C + 2 9 2/3 v 2/3 ( + 2 9 

Since the initial condition is n ~ 5(z/ — i/q) one can suppose that 

-^-v 2/3 n « v 2/3 -^-n 
av av 

Then the kinetic equation (J2J) can be written as 



an 2/3 a c + 2 a 2 _ 1/3 c + 2 d 

— — = —v 1 -^—{Cn -z-n) v • [in — 

dt <9z/ VS 2 dv ' 3 VS 2 

The second term 

II = A u -^ {(:n _(±ld 

3 vs 2 dv 1 

can be comparable with the fist term 



n 



only when 



, 9 C + 2 d .. C + 2 3 

— - — — Ti) — f — — n\ 

'<9z/ S 2 dv Jl IS 2 <9z/ 1 

which means that we are near the stationary solution corresponding to 

(Cn - ~l 2 ^ n) = J s = const 
2 dv 
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where J s is the stationary flow. But here we consider the growth of a single 
droplet, the situation is opposite and the true initial condition is 

n(v ) ~ S(v - v ) 

So, here the second term can be neglected and the relation Q can be justified. 

Now we can see that having neglected // in kinetic equation one can 
reduce it to 

dt 3dr [ ^ 2 3r 2 «9r } K ) 

for r = z/ 1 / 3 . 

We see that the last equation is the ordinary diffusion equation with 
constant velocity of regular growth and the coefficient D r of diffusion along 
r-axis written as 

°r = ( 5 ) 

Since D r is a decreasing function, one can see that the distribution function 
over v (as a function of r) will be a well localized function in r-scale (and, 
hence, in z/-scale). So, the characteristic relative width of n will be small. 

This conclusion will also lead to the self consistency of negligible character 
of the second term in kinetic equation. 

Under the constant value of D r the characteristic width 5r will have an 
order r 1//2 and the characteristic value of r is many times greater than the 
characteristic width. Here the situation leads to a more strong inequality. 

To see the effects of stochastic growth one can also calculate the derivative 
< ±v >. The expression is the following 

d d 



Then 



d d r°° 

< —v >= — / vn{v,t)dv 
dt dtJo K ' 

d r,d 2/3 C + 28n. 

or 

d f°° >. 2/s , ( + 2 f°° , n/ndn 

< — v >= - / n(v 2/3 dv + / cW 2/3 — 

at Jo 2 Jo ov 

Note that we have to use namely precise equation (J2J) without transformation 

(0J). If we use ® we get the additional wrong factor 5/3. 

The first integral 

h = / nCv^dv 
Jo 
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corresponds to the regular growth and the second integral 

r C + 2 /-°° , /3 <9n 



/■OO 

/ cW 2 / 3 
Jo 



<9z/ 



represents corrections. 

One can see that integration by parts gives 

2C + 2 {<*> _ 1/3 

Since n is well localized in a relatively small region I 2 can be presented as 

J 2 ~-_<„ > 
and it is rather small in comparison with I\ which can be estimated as 

h ~ c < ^ 2/3 > 

So, the smallness of correction to the mean coordinate is proven analyti- 
cally. 

One has to note that the mentioned conclusion is rather typical for the 
processes of such kind. We can make our arguments more strong if we recall 
the Uhlenbeck process. One can note that for the Uhlenbeck process 

dn d . . d 2 
_ = 7 _( a n)+iJ_ n 

with two parameters 7 and D the Green's function is known and has the 
exponent form: 

nix t\x> t>) - ! 7 exp( 7(x-exp(-2 7 (t-Q^ 

1 ' 1 ' J " V 2wD(l - exp(-2 7 (t - t'))) P[ 2D(1 - exp(-2 7 (t - t'))) } 

This exponent form lies in correspondence with the absence of correction to 
the regular growth. Really 

since dx/dt ~ 7 x. Then 

J xn(x)dx = Xq + J \x — x )n(x)dx 
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with arbitrary x . 

The coordinate of the regular growth xq will be exp(— 27 (£ — t')x'). Then 
the Green's function is the Gaussian symmetrical on x— xq. Then < x >= xq. 
As the result the first term x = dxo/dt describes the regular evolution and 
the second term J(x — x )n(x)dx is absent. This is one of the reasons why 
the Green's function is known here. 

In our process the analogous property is absent. But the last example 
shows that at least the effects of the difference between d < x > jdt and 
< dx/dt > are moderate because in situations dx/dt ~ const and dx/dt ~ x 
(our situation lies between these laws) they are absent. 

4 Numerical results for dispersion as a func- 
tion of the final size 

The next task is to calculate dispersion Su as a function of Uf in and u . Here 
there is no need to average over such a giant number of attempts, it is quite 
sufficient to take 200 attempts to get suitable result for dispersion. The 
results of simulation are the following: 

• For u = 20 one can see in the Figure 1 the following dependence of Su 
on u fin 

One can see that it can be approximately the straight line 

Su = 42 + 0.095u fin 

• For v>o = 50 one can see the analogous dependence of Su on Uf in drawn 
in Figure 2 

One can see that it is also approximately the straight line 

Su = 17 + 0.094z/ /in 

The result of these two numerical pictures is that one can not say that 
the fluctuations disappear when u goes to infinity. The result is that 

Urriv-HxjSv = 0.095z/ 

The direct sequence is that one has to take fluctuation effects of growth 
into account. 
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Figure 1: Dispersions for u = 20 as a function of x = Uf in 

Fortunately, the true result is a little bit more optimistic. In reality bvjv 
becomes small for enough big v. More precisely the result of mentioned 
simulations give 

5u ~ v l ~ e 

where 

e ~ 0.22 ±0.05 

It can be seen from corresponding pictures drawn in logarithmic scales. 
The following Figure 3 shows the dependence of In 5u on v for u = 20 
The linear approximation 

ln<5z/ = 0.741nz/-0.76 

is also drawn. 

The same dependence for uq = 50 is drawn in Figure 4 
The best linear approximation 

]n6u = 0.83 In i/- 1.03 

is also drawn. 

So, the convergence bvjv — > formally takes place. But for every concrete 
(finite) situation one has to take this correction into account. 
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Figure 2: Dispersions for u = 50 as a function of x = Ufi 



5 Theoretical explanation for dependence of 
dispersion on the final size 

Here we are interested in the dependence of 5v on v^ n . This dependence will 
be like a power one and the value of the power has to be determined. 

The starting point of explanation will be equation (jlj). We can see that 
D T = 3^ has a singularity at r = 0. Certainly, r = is out of range of 
consideration. At least D r is a rapidly decreasing function. It means that 
the main blurring takes place in initial period of time of growth (i.e. some 
$initiai r has been formed) and later this blurring is mainly translated along 
z/-scale. 

It is rather easy to estimate the role of translation. Really, due to trans- 
lation 

8v = Sinitiair (— ) = 5 initia ir 3z/ 2/3 



dr ' 



since the value of initial dispersion 5r is fixed we come to the power like 
dependence 



X 2 /3 
™ ~ "fin 



The power is 2/3. 



12 



5 6 7 8 9 10 

x 



Figure 3: Logarithm of dispersions for u = 20 as a function of x = lnz/^. 



One has also to consider the standard diffusion of purely independent 
events with blurring in the r scale. This leads to 5r ~ r 1 / 2 . Then 

5u ~ r 1 ' 2 = u 1 ^ 

The power like dependence together with scaling invariancy means that 
some fixed part of evolution forms the dispersion 5i n iu a ir. For this value one 
can take 

c 1/2 _ 1/6 

initial 11 ~ r initial ~ ^initial 

Since initial = ocUf in , where a is some fixed parameter, one can come to 

^initial" ~ ^fin 

and 

1/6, dv 2/3+1/6 0.83 

Numerical simulations confirm this result. 
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Figure 4: Logarithm of dispersions for v Q = 50 as a function of x = \nuf in 

6 Numerical results for dispersion as a func- 
tion of the initial size 

The problem of establishing the dependence of 5v is not yet solved because 
there exists the uncomfortable dependence of Su on the starting size vq. 

One can not put u to zero and solve this problem because in numerical 
simulation one has to exclude the dissolution of embryos. 

The simulations with finite is shows the irregular behavior of 5v for small 
v. Seeking for the power like dependence we shall draw curves in logarithmic 
scales. The dependence of In 5u on In v is drawn in Figure 5 for uj in = 1000 
and Figure 6 for i>f in = 2000 

The straight lines drawn in these figures are the following 

• In Figure 5 

y = 6 — 0.3a; 

• In Figure 6 

y = 6.3 - 0.27x 
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Figure 5: Logarithm of dispersions for h>f in = 1000 as a function of x = In u 

7 Theoretical explanation of dependence on 
initial size 

Again we are seeking for the power in the power like dependence of 5v on vq. 
Having noticed that D r ~ z/~ 2 / 3 we can rewrite this equation as 



where t' is the moment corresponding to r. 

We see that D r is mainly concentrated in initial moments of time. This 
allows us to state that approximately 



Having calculated the last integral and interested in behavior at the lower 
limit of integration we come to 



D r ~(t- t') 



-2 




r 1 / 2 ~ 5r 



In terms of v ~ t 3 



one can rewrite the last equation as 
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Figure 6: Logarithm of dispersions for Uf in = 2000 as a function of x = In u 

The numerical results confirm in general features the small negative value 
of the power. But the some discrepancy between theoretical model and 
results of simulation still remains. 

Then we have to use the ideas of the self-similarity. There is no difference 
between v and Uf in . Really, the value z/ is the final value for the region 
[0, uo] which is out of consideration. The dispersion was supposed to be like 
v\j 2 . But since u can be interpreted as a final size, we know that dispersion 

1+2 

really will be z/ 6 3 . 

Then we have to increase the power of the previous result in 

(l/6 + l/3)/(l/2) 

times. This leads to 

Sr ~ i/-[( 2 /3+V6)/(l/2)]*(l/6) = ^-5/18 

Then the coincidence between theoretical and experimental results be- 
comes satisfactory. 
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8 Calculations for the kernel 



Since the functional dependence on vq and on Ufi n have been determined we 
don't know only the numerical coefficient in the formula 

r , -5/18 5/6 

bv = constu v/ in 

This coefficient can be determined by the unique computer simulation or by 
application of analytical model to one of the concrete situations. 

Here we shall show no more than an illustrative example of calculations. 
The weak feature is that it does not correspond to all approaches used above. 

Let us choose the situation used in simulations for the mean value. The 
value of dispersion has to be calculated as following: 

• The moments of time 

to = 3^o V3 /C (6) 



(7) 



and 

t 

have to be determined. 
Dispersion in r scale is 

5r = v/2N total 

where N tota i is the effective total number of events. It cam be approxi- 
mately calculated according to the following formula 

Ntotai = 2 I' 1 D r dt 

J to 



to have an analogy with br = yADt in a case D = const. Having used 
© and keeping in mind 5r ~ y/ADt we come to 



V Jto 9 2 C 2 t 2 V vs >C t t s ' 



Dispersion in v scale will be 

5 

After substitution of © and (JJJ) we come to 



Su = 5r * 3i$J 



bv = 3u 



2/3 
fin a 



C + 2, C 



\ C 2 v 3z. 1/3 3u 



fin 
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• It is quite evident that the dispersion of self blurring 

S u = p(C + 2)/Cf/m 

has not to be added. This is because we have not divided the whole 
region into two parts but integrate D r over the whole time interval. 

The result for dispersion without self blurring will be 

5v = 155 



The addition of dispersion of self blurring leads to 

8v ~ 200 



It practically coincides with the result of simulations. 
As the result the dependence of 5u will be 

Su = l55(^) 5 / 6 (-) 5 / 18 
V 1000 ; V zV 



9 Application for concrete systems 

Decay 

Now we shall apply the already known formulas for concrete systems. It 
is known that two typical types of external conditions are the decay of the 
metastable state [Sj and the so-called dynamic conditions for condensation 

Hj. 

The global theories of evolution based on averaged rates of nucleation 

(formation of droplets) and droplets growth were given in [3] for decay and 

in [3] for dynamic conditions. 

The peculiarities of stochastic appearance of droplets were investigated in 
|S] for decay and in [Hj, [Zj for dynamic conditions. In these publications 

the rate of growth for supercritical embryos of liquid phase is supposed to be 

the regular one. 

The results in description of stochastic effects in [2], [S], 0, [I] were 
achieved by application of algebraic structure of the nucleation period. This 
structure will be very important below. 

In the situation of decay the structure was the following one: 
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• Until the moment when the prescribed (by recipe of the monodisperse 
approximation) part of the total (averaged total) number of droplets 
has been appeared the system simply waits. This is the first subperiod 
of the nucleation period. 

• Then the second subperiod begins. In this period the rest droplets 
appear. The vapor consumption and, thus, the behavior of supersatu- 
ration is fully governed by the growth of droplets appeared in the first 
subperiod. 

• Dispersion of the number of droplets appeared in the first subperiod 
equals to zero, the system is simply waiting for the appearance of the 
necessary amount of droplets. 

• Dispersion of the rest droplets, i.e. of the number of droplets appeared 
in the second subperiod equals to the standard dispersion. But due 
to the fact that the number of the rest droplets is less than the total 
number of droplets the final dispersion will be less than the standard 
one. 

The number of droplets appeared in the first subperiod will be marked 
by N first- The number of droplets appeared in the second subperiod will be 
marked by N rest . 

All conclusions except the last one remain valid also under the stochastic 
model for the droplets growth. The last conclusion has to be reconsidered. 
Really, earlier the cut off by the vapor consumption was absolutely regular 
since the number of droplets in the first subperiod has to be absolutely fixed 
and, hence, the action of these droplets was absolutely prescribed. 

Now the situation is changed. The cut off of the rest part of spectrum 
can be initiated in different moments of time since the droplets from the first 
part of spectrum (let us call them as the the leading droplets) can grow faster 
or slower. 

One can illustrate the effect of the different velocity of growth by the 
following picture 

• Suppose that the droplets grow slower than the averaged velocity. Then 
the moment of the cut off will be later and more droplets than the 
standard value appear in the rest part of spectrum. 
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• Suppose that the droplets grow faster than the averaged velocity. Then 
the moment of the cut off will be earlier and less droplets than the 
standard value appear in the rest part of spectrum. 

We shall call this picture as the "unbalanced effects of growth" (UEG) 
From the effective monodisperse approximation [5] we know the "initial" 

size and the final " size" . The number of droplets in the monodisperse peak 

in also known. 

Certainly, the picture of UEG isn't absolutely correct. It is clear that 
the balancing factor appear. This factor is similar to that which leads to the 
absence of dispersion in formation of the leading droplets. 

Suppose that the initially determined leading droplets grow slower than 
ordinary. Does it mean that only these droplets have to be taken into account 
in vapor consumption during the nucleation period? Certainly not. We 
have simply to include into the group of leading droplets some new droplets 
appeared a few moments later than the initial boundary between the leading 
droplets and the rest droplets. So, the existence of a balancing force is evident 
here. 

Now suppose that the initially determined leading droplets grow faster 
than ordinary. Then we have to exclude from the group of leading droplets 
some the last droplets appeared just before the initial boundary between the 
leading droplets and the rest droplets. So, the balancing force here also takes 
place. 

It is clear that the picture with unbalanced effects of growth will estimate 
the stochastic effects of growth from above. We shall use this property to 
calculate the stochastic effects in this picture and then to see that they are 
small enough. 

In reality the result from stochastic growth will be some times smaller. 

Now we turn to the general case of asymptotic of the big number of 
molecules in critical embryo. 

To grasp the stochastic effects from above we use the picture of unbal- 
anced effects of growth and calculate the effects from the stochastic appear- 
ance of droplets and from stochastic growth into dispersion. We suppose that 
one can add dispersions caused by both stochastic causes. It is true at least 
when one of dispersions is much smaller than another. Namely this situation 
will take place. 

Dispersion of appearance is 
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and dispersion of growth will be connected with the fact that the number of 
the rest droplets N rest will fluctuate. 

If we forget that N rest is the argument of the square root which smoothes 
the result we will get the estimate from above. So, we will forget. 

So, we will calculate the dispersion 5\ of one droplet caused by stochastic 
growth from the size u st up to the final size. 

Fortunately the initial size is well determined in the advanced monodis- 
perse approximation and it is not zero. In the previous initial monodisperse 
approximation it was zero and the consideration of stochastic effects required 
to give in j^j a new version where the initial position of monodisperse spec- 
trum (the z -coordinate) in renormalized scale is 0.4 — 0.336. Namely the 
new version has to be used. 

So, the value of Vq for all droplets in the monodisperse spectrum will be 

z/ 1/3 = (0.4 - 0.336)^ 

One can also determine in a following simple manner. The actual 
height of monodisperse spectrum is the unperturbed rate of nucleation (of 
appearance) Jo- Since 

N fi rs t ~ J 5r 

we can get the characteristic width of the spectrum 5r. This value will be 

1 /3 

initial value for v . It is clear that approximately this way leads to the 
same value of vq. Then since the dependence on vq is like v Q 1//6 we see that 
the results of two ways to determine v$ will be approximately equivalent. 

The zero value of vq would lead to divergence, now such a danger is over 
and one can see that the asymptotic at v — > oo for characteristic relative 
deviation in v is 

r 5/6 / 5/18 

OV = » fin/ » fin 

and goes to zero. The total characteristic relative deviation will be 



for one droplet in the monodisperse peak 

n 5 /6 / 5/1 

L>unbal ~ V 'fin I v 'fin 

for all droplets in the monodisperse peak 
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It will be many times less than the relative dispersion 



stoch 




of the number of droplets caused by stochastic appearance. 

We have to note the absence of the shift of the averaged value of the final 
size from the value calculated on the base of the regular law of growth. The 
order of this deviation is less than the characteristic attained value. If this 
shift was essential the action on the averaged number of droplets would be 
significant. 

One has to note that the in last remark the regular law means the regular 
discrete law and the main effects will be the effects of discrete character of 
the droplets growth. 

Dynamic conditions 

In dynamic conditions all previous arguments are valid. The initial size 
vq has to determined by the following way. Since in the approximation of 
monodisperse consumption [7| the monodisperse spectrum is formed already 
at z — — 3/c one can say that it is formed at the ideal conditions (the real 
supersaturation coincides with the ideal one). Then the characteristic width 
of spectrum will be c _1 . Then the starting value value u^ 3 will have the 
order c _1 

1/3 -1 
Vq ~ c 

The final value u 1 / 3 will have the order (3 -v- 4)c _1 



All other constructions are absolutely analogous. 

Certainly, the numerical simulation can not give correction terms where 
the effects of stochastic growth can be noticed because these effects will be 
corrections for corrections to the main result. The corrections connected 
with the discrete growth by steps of absorption of one molecule will be more 
significant. 

One has also to note that at the further periods of evolution (after the 
end of nucleation) the stochastic effects of growth will be more important. 
This consideration will be presented in a separate publication. 



v 



,V3 
fin 



(3-4)c 



,-i 



22 



References 



[1] Kurasov V., The form of embryos sizes in the first order phase tran- 
sitions, Theoretical and Mathematical Physics, vol. 131, number 3 , 
p.503-528 

[2] Kurasov V., Preprint cond-mat/0410616, Explicit two cycle model in 
investigation of stochastic effects in diffusion regime of metastable phase 
decay, 61 p. 

[3] Kuni F.M., Novojilova T.Yu., Terentiev I.A., Kinetics of heterogeneous 
nucleation under the instantaneous creation of initial supersaturation. 
Russian Journal of Theoretical and Mathematical Physics, v. 60, p. 276, 
(1990) 

[4] Kurasov V.B., Phys. Rev. E, vol.49, p.3948 (1994) 

[5] Kurasov V., Preprint |cond-m at/04l 6~774| Monodisperse approximation 
as a tool to determine stochastic effects in decay of metastable phase, 
19p. 

[6] Kuirasov V. , Preprint cond-mat / 0412142 , Corrections to a mean number 
of droplets in nucleation, 27 p. 

[7] Kurasov V., Preprint cond-mat/0412141, Dispersion of nucleation under 
the smooth variation of external conditions, 25 p. 



23 



